Adiabatic information transport in the presence of decoherence 



I. Kamleitner, J. Cresser, and J. Twamley 
Centre for Quantum Computer Technology, Physics Department, 
Macquarie University, Sydney, New South Wales 2109, Australia 

We study adiabatic population transfer between discrete positions. Being closely related to 
STIRAP in optical systems, this transport is coherent and robust against variations of experi- 
mental parameters. Thanks to these properties the scheme is a promising candidate for transport 
of quantum information in quantum computing. We study the effects of spatially registered noise 
sources on the quantum transport and in particular model Markovian decoherence via non-local 
coupling to nearby quantum point contacts which serve as information readouts. We find that the 
rate of decoherence experienced by a spatial superposition initially grows with spatial separation 
but surprisingly then plateaus. In addition we include non-Markovian effects due to couplings to 
nearby two level systems and we find that although the population transport exhibits robustness in 
the presence of both types of noise sources, the transport of a spatial superposition exhibits severe 
fragility. 
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I. INTRODUCTION 

The coherent transport of quantum information is an 
essential clement in any scalable architecture for a quan- 
tum information processor. The interconversion of static 
qubits into flying qubits is difficult and may not be a 
feasible solution for the coherent spatial transport of 
quantum information. As a possible alternative much 
attention has been focused on dark state adiabatic pas- 
sage for coherent state quantum transport. Originally 
studied in the context of quantum optics [2] where it is 
called stimulated Raman adiabatic passage (STIRAP), 
dark state transport uses the existence of a "dark state" 
which is a zero-energy cigcnstatc of a driven quantum 
system. By manipulating the driving of the system one 
can sculpt this dark state to coherently transport quan- 
tum states using STIRAP-like procedures. This intra- 
atomic dark-state transport has been demonstrated ex- 
perimentally [3 . However the method has more recently 
been applied to spatial transport of quantum information 
(which we specifically denote CTAP - coherent transport 
by adiabatic passage following j4j), in a variety of physical 
systems including chains of neutral atoms [5], quantum 
dots [HEBE], superconductors [8 j; and photons in nearby 
waveguides [§]■ It has also been proposed as a crucial 
element in the scale up to large quantum processors |10j . 
The method possesses two very crucial benefits over other 
quantum transport methods: since the transport is via a 
zero energy state the quantum state acquires no dynam- 
ical phase and due to the adiabatic theorem, the process 
is very robust to a wide range of system variations. How- 
ever, an important question, particularly with regard to 
the use of CTAP within large scale quantum computer 
architectures, is to determine the effects of decoherence 
on the transport. The effects of dephasing and sponta- 
neous emission has previously been examined in the case 
of STIRAP dark-state transport in a three level atom in a 
A configuration [HI [12] . In that work a master equation 
was postulated and its effects on the transport studied. 



In our work we examine the effects of two types of 
physically-motivated decoherence sources effecting the 
CTAP transport in a quantum chain. We first study 
the effects of delocalised measurements on the systems 
making up the CTAP chain. In particular we imagine 
an electron on a chain of quantum dots (QD) and each 
is measured by a nearby quantum point contact (QPC). 
These QPCs however, are non-local measurement devices 
in that their charge sensitivity falls off with distance. 
Such devices or similar will be required to cither mod- 
ulate or readout the quantum information in each quan- 
tum dot in the CTAP chain in a real device. In large 
scale quantum processors one will routinely wish to have 
quantum information in a superposition of two "distant" 
spatial locations. It is known that from numerous stud- 
ies of cat-states in quantum Brownian motion - a single 
harmonic oscillator coupled to a bath of harmonic oscil- 
lators - the rate of decoherence suffered by the cat grows 
quadratically with the spatial separation of the two su- 
perposition states of the "cat". We find that such an 
effect is also present in our case, i.e. the decoherence 
rate of a "cat-state" on the CTAP chain increases with 
cat-seperation, but surprisingly we find that this decoher- 
ence rate saturates beyond a critical cat spatial sepera- 
tion. This is a positive result for the CTAP transport 
protocol and is essentially due to the rapid spatial fall off 
of the measurement sensitivity of the QPCs. 

This first model is an example of Markovian deco- 
herence. We also include a second non-Markovian de- 
phasing source, and consider that each quantum dot 
interacts with a nearby two level system (TLS). This 
could model unlocated two-level fluctuators in a solid 
state CTAP scheme. Interestingly we find that qubit 
transport still seems relatively robust in the presence of 
these combined Markovian and non-Markovian decoher- 
ence sources. More worryingly however we find that this 
non-Markovian dephasing slightly entangles the two level 
systems with the transported qubit. This effect does not 
seriously detract from the transport of an electron iso- 
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FIG. 1: Realization of non-local measurement as proposed in 
|16| . The electric current through a QPC is influenced by the 
appearance of an electronic charge in its near environment. 

lated on one QD but causes serious degradation if the 
qubit transported is in a superposition state. As the 
latter situation will be the typical case in a large scale 
quantum processor the present analysis might indicate a 
much lower density of unlocated TLSs will be required 
when using CTAP in large scale quantum computer. 

We proceed in section [IT] to motivate and model two 
sources of decoherence and solve the corresponding mas- 
ter equation without interactions between different QDs. 
Transport of the electron by CTAP is described in sec- 
tion IIIII under both the Markovian and non-Markovian 
noise models. Interestingly we are able to analytically ob- 
tain an equation describing the fidelity of the Markovian 
dynamics while we numerically solve the non-Markovian 
case. It is here where differences between Markovian and 
non-Markovian dephasing become most apparent. We 
summarize our results and discuss their relevance in sec- 
tion us 

II. MODEL 

A. Measurements 

As we mention above, we first consider each of the 
CTAP quantum systems to be continually measured. 
This gives rise to Markovian decoherence. We take the 
specific example of an electron on a chain of quantum 
dots and consider the measurements to be made by quan- 
tum point contacts (QPC) situated close to the quantum 
dots as shown in Fig[l] The measurement executed by a 
QPC is caused by the modulation of the current trans- 
mitted through a QPC due to the presence of a nearby 
electron. The QPC current is modulated by a factor 
1 — f , where r is the distance between electron and QPC 
and a is a constant reflecting the properties of the QPCs 
(see [IB]). This modulation results in an indirect position 
measurement of the electron's spatial position on the rail 
of QDs. However, it is a non-local measurement because 
even an electron on the neighboring QD influences the 
transmission through a QPC. The localness is parame- 
terized by 2> i- e - the distance between QPC rail and 
QD rail over the distance between two neighboring QDs, 
with small values representing more local measurements. 
Furthermore — parameterizes the sensitivity (signal over 
noise) of the measurements and is typically small and 
hence the measurements are weak ones. Such measure- 
ments are properly described in the language of positive 



operator valued measurements (POVM) [TT1 US] , 

The purpose of a measurement apparatus is the read- 
out of quantum information for which one would like to 
use strong local measurements. This can be approxi- 
mated by using a large number of weak, non-local mea- 
surements of the type described. A large number of mea- 
surements in a reasonably short time is achieved by using 
a high measurement rate which in turn can be realized, 
for instance, by applying a voltage to the QPCs. How- 
ever, such measurements also act as a source of decoher- 
ence. During quantum unitary operations such as trans- 
portation by CTAP it is preferred that this decoherence 
is absent. However this might not be completely achiev- 
able in practice and one may be left with a small deco- 
herence rate due to the coupling to the QPCs during any 
quantum operation. Thus measurements as a source of 
decoherence are included in the analysis presented here. 

We restrict our treatment to the case of having only 
one electron in the rail of QDs. We also assume that the 
electron can only occupy the ground state of the QDs, 
and we take \i) A , to be the quantum state of the electron 
in the i th QD. Furthermore we neglect all interactions 
depending on the spin of the electron. Then, {\i) A , i — 
1, ..,iV}, form a basis for the Hilbert space Ha of this 
electron on the QD rail with N dots. In the following 
we take the limit of a long rail, N — * oo. We denote the 
distance between the i-th QPC and the j-th QD by 

rn = v/« 2 + (N-jM) 2 (i) 

where a and d are defined in Fig [T] 

The probability of the j th QPC detecting the presence 
of an electron on the QD rail can be written as [T7] 

P j (p A )=Trfap A ) (2) 

where p A is the state of the electron on the QD rail and 
TTj is the effect operator corresponding to the QPC mea- 
surement at site j. If the electron is spatially localised 
to be only on the i th QD, i.e. in the state \i) A , Eq. ^ 
reduces to 

Pji\i) A (i\)= A (i\*j\i) A - ( 3 ) 

As we noted above the measurement sensitivity of the 
QPC depends on the distance r%j. The presence of an 
electron a distance r,j away from the QPC decreases the 
current flowing through the QPC by a factor 1 — — and 

this leads to a reduced detection probability, 

PiQi) A (i\)<xl-£- . (4) 

Fulfillment of Eq. @ is certainly achieved with the mea- 
surement effect operators 

N / \ 

i—l x J / 



3 



with 



N 



N- 



1 



N 

Ea 
Nra 



(6) 



to satisfy Ylj=i n j ~ Note that each effect operator 
is almost proportional to the unit operator which reflects 
the weakness of the measurements being performed. 

Measurement theory states that the transformation of 
the density operator due to a measurement is described 

by PUCES] 



A 3 P A A] 



(7) 



with Aj — Uj^/Wj for some arbitrary unitary operators 
Uj . This unitary depends on the interaction of the quan- 
tum system and measurement apparatus during the mea- 
surement process. In our case, if the electron is at site i, 
its position will not be changed by a detection event by 
a QPC at site j. For simplicity we also assume that the 
measurement does not introduce relative phases within 
the rail of QDs giving Uj = 1 and A — A\ which means 
the measurements influence the electron state as little as 
possible. 

To derive the master equation we now assume that 
detection events in the QPCs occur uniformly at random 
and at a constant rate R. Following [H5], we can write 
down the master equation describing the evolution of the 
density operator as 



dp A 
dt 



N 



i[H A , p A ] - Rp A + rJ2A jP a A j 



(8) 



This equation is in Lindblad form with Lindblad oper- 
ators Aj and thus the evolution is Markovian. For the 
case when Ha involves no electron transport along the 
QD rail, Eqn. Eq. (|8) possesses a stationary state which 
is diagonal in \i) A , and thus the evolution corresponds to 
a pure dephasing type of decoherence. We note that R 
should scale proportional to N as each QPC contributes 
equally to the overall measurement rate. 

For now we take Ha = 0, but later we introduce a 
time dependent Hamiltonian to induce CTAP. Express- 
ing Eq. Q in the basis \k) A , we obtain p A k — 0, for 
diagonal entries and p A , — —T^ip^i for off-diagonal ones. 
The decoherence rate 
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(r^ is defined in Eq. pi) can be shown to saturate for 



large distances \k — l\cT^> a. To simplify these formu- 
las one can use the limit of non-sensitive measurements 
^- < 1 which should be well justified in experiments [15] . 
This weak measurement limit corresponds with the in- 
ability of the measurements to give detailed information 
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FIG. 2: As in Fig [T] but with additional local coupling to 
two-level systems. 



on the position of the electron in the QD rail. Weak 
measurements of this nature feature in many models of 
continuous monitoring of a quantum particle's position 
20 . In the limit of weak measurements we obtain 



D 2 N 

Ra z x - 
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(10) 



which, for spatial separations larger than the threshold, 
\k — l\d 3> a, limits to 



Ra 2 coth (na/d) 
4Nd 2 a~Jd ' 



(11) 



We also note that we can execute local measurements, 
which give maximal information about the particle's po- 
sition in the limit ^- = Sij, while keeping a finite, e.g. 

in Eq. In this case we find 



R 



T kl ='-(2 + a- 2VY+a) 



(12) 



The saturation of the decoherence rate in Eq. (Ill, is 



somewhat surprising when one compares this with the 
similar situation for a free particle in a spatial superposi- 
tion cat-state, experiencing continuous position measure- 
ments |21j . In that case the decoherence rate suffered 
by the particle increases without bound according to the 
spatial separation of the cat, i.e. T XltX2 oc (xi — X2) 2 . 



B. Coupling to TLSs 

We now consider a further source of decoherence. It 
is highly likely that in any physical device there will be 
unknown accidental two level fluctuators nearby to the 
quantum dot rail (see Fig. |2|. If these unknown two 
level systems (TLS), can couple to the electron on the 
rail then these systems effect a source of decoherence 
which exhibits memory, i.e. is non-Markovian. Quan- 
tum coherences on the quantum dot rail can be trans- 
ferred to the nearby TLSs, where they can remain for 
a period, before being transferred back. Typically the 
analysis of these types of non-Markovian effects are com- 
plex but in the following we are able to derive analytic 
solutions of the resulting reduced dynamics of the quan- 
tum dot rail. For simplicity we assume these fluctuators 
have no internal dynamics other than their coupling to 
the quantum dot rail [23] . Furthermore we assume that 
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these TLSs do not experience significant decoherence on 
the time scales of the transport. We are aware that these 
assumptions may not be completely satisfied in all realis- 
tic situations. However this model will serve to highlight 
the striking difference between the Markovian and non- 
Markovian evolutions in, for instance, the transport of a 
spatial superposition. 

We use | IV- and |0) ■ as basis of the Hilbert space Hj of 
the j—th TLS such that the interaction Hamiltonian is 
diagonal. If the electron is on the j-th QD, it is assumed 
to induce a phase shift on the j—th TLS, so that the 
Hamiltonian acting in the product Hilbert space H = 



Ha <8) 7 -=i H-j °f electron and TLSs reads 



N 



Hint — ^ 



Xn \n) A (n\ <8> 



<r*,n Qs) h 



(13) 



The coupling constants Xj are considered to be constant 
in time, and a z j and lj are the Pauli Z-matrix and the 
identity operator, respectively, acting in Hj. 

As typically assumed, we now take the initial state to 
be in product form p(t = 0) = p A (8 Pq LSs , where Pq LSs 
does not have to be a product state of the different TLSs. 
We can now express the master equation describing the 
time evolution of the density matrix of the combined QD 
rail and TLSs under the effects of the above measure- 
ments (see Eq. ((8|), to be 



dp 
dt 



-i[H int ,p] - Rp 



N 
3 = 1 



(14) 



with Aj = 7tF-(8> 1® N - 

After some effort, one can trace out the TLSs to find 
the non-Markovian master equation for the reduced den- 
sity matrix of the QD rail for an arbitrary initial product 
state of the TLSs, as 



dp A 
~dt 



N 

-R P A + (RA n p A A n - iA n {t) [\n) A {n\y 



n=l 



+ Jn (t)[\n) A (n\,[p A ,\n) A (n\}}), 



(15) 



where the last two terms describe the effects of the TLSs 
on the QD rail. Here the definition 

j n {t) - iA n (t) = Xn — . 7 (16) 

COS Xnt + «W„ Sin Xnt 

is used and Lo n =Tr[p„<7 z .„], is the inversion of the n-th 
TLS. Decoherence due to the TLSs is described by j n 
whereas A„ represents the Lamb-shift. If u n — ±1, then 
A„ = ±Xn and j n — 0. In this case the decoherence due 
to the coupling to TLSs vanishes and the coupling to 
the TLSs results only in a change of the energy of states 
\ n ) A by i%n- However, a more interesting case is when 
u> n = 0, since it includes the case where the TLSs may 
initially be in a complete mixture /O„(0) = \ 1. In this 



case we find A„ = 0, and 7„ = Xn tanx n t. The presence 
of singularities in j n (t), may cause difficulties in studying 
the time evolution of p A , using normal methods. To avoid 
this we do not work directly with ( 15 ). Instead we return 



to solve the complete dynamics of the coupled QDs and 
TLSs and then trace out the latter to obtain p A {t). 

To achieve this we first we go to the interaction picture 
by transforming Eq. (14 1 with 



N 

e *H int t = |n)^(n|(g)(cos(x n t)l„ + ^sin(x„^)cT ^ ,„) (g) l r 

n— 1 j¥" n 

(17) 

The combined density operator in this picture is given by 
p(t) — e lHintt p{t)e~ lHintt , and is governed by the master 
equation 



dp 
dt 



N 



= -Rp + R^2 AjpAj . 



(18) 



Note that e tHi "'* Aje^ 11 ^ = A r Taking the compo- 
nents of Eq. (18 1 in the QD Hilbert space, /5 fci (t) = 
A.{k\ ~p \l) A (which is still an operator in the TLSs Hilbert 
space) , we find the solution is given in a similar manner 
to that of the previous subsection to be 



.-Tut 



p M (0) = e-^puiO) 



(19) 



where T^i is given in Eq. (j9j. After transforming back 
to the Schrodinger picture and tracing over the TLSs we 
find for the components of the reduced density operator 



Pkki*) 



-T M t 



cos (xkt) cos (xit)Pki(0) 



(20) 
(21) 



where we have assumed that the initial states of the TLSs 
is a completely mixed state 



Pi(0) 



1 



1. 



(22) 



Hence we see that information lost to the TLSs can return 



to the rail via the oscillatory terms in Eq. (21 1, which 



is in contrast to Markovian decoherence induced by the 
measurements. Note that the non-Markovian behavior 



of Eq. (21 1 results from tracing out the environmental 
TLSs. 



III. CTAP 

Coherent Tunneling by Adiabatic Passage relies on the 
adiabatic theorem in the sense that the electron will al- 
ways be in an eigenstate of the time dependent Hamilto- 
nian. For the transport from the m-th to the n-th QD, 
the Hamiltonian is designed such that there exists an 
eigenstate which is \m) A at the beginning of the trans- 
port, and then changes continuously to |n) A during trans- 
port. This way of population transfer has the important 
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property of being robust against small variations of the 
Hamiltonian and the transport time. This is crucial in 
many experiments since these parameters are often hard 
to control. The drawback of CTAP is a relatively long 
transport time which is limited by the adiabatic theo- 
rem and usually is about an order of magnitude longer 
compared to diabatic population transfer. 

To transport an electron from position m to position n 
we assume that we can control tunneling rates between 
neighbouring QDs, but we do not require any next to 
neighbour coupling. The Hamiltonian then reads 



H A = 



( o n m o • 

o n m+1 o 



ooo 
ooo 
V o o o 



o o o \ 
ooo 





o n„- 2 o 

ft n -2 Q n _i 

■ fin-! / 



(23) 



The time dependent manipulation of the tunneling 
rates flj can be achieved in an experimental setup via 
external gates as is outlined in |J. The tunneling rates 
Cl m = Qp and fl n -i = are often called pump and 
Stokes pulse, respectively. The application of CTAP can 
roughly be outlined as follows: 



1. While keeping fl 
j = m+ 1, . . . , n- 



= one switches flj 
1 to ensure all energy eigenstates 
for an electron situated between m and n are non- 
degenerate. 



2. Then £l m is increased and fi n _i is decreased until 
it vanishes. This process has to be done slowly to 
satisfy adiabaticity. In this step the electron moves 
from the m— th to the n— th QD. 

3. Finally all couplings can be set to zero. 

To avoid any geometric phases one can use real and pos- 
itive tunneling rates flj. In the following we will refer 
to these steps as step one, step two and step three. We 
emphasize that the electron moves exclusively in the sec- 
ond step and therefore the first and the third steps can 
be done arbitrarily fast (up to experimental limitations) . 
Hence we set t = to at the beginning of step two and 
t = t max at the end of step two. It will be assumed that 
at t = to the probability of finding the electron on the 
j-th QD with j = m + 1, • • • , n is zero. This is reason- 
able since the adiabatic theorem does not hold for these 
states and we want to transport the electron from the Tro- 
th QD. Furthermore we use that n — m is even, because 
only then all states are non-degenerate after step one. 



A. Measurements 

To introduce our technique for solving the master equa- 
tion we first neglect the coupling to the TLSs. At 



1) eigenstates of Ha which are su- 



/ A ' 



are denoted by \ipj(to)) 



t = to the (n — m 
per positions of \m) 
with j = ^f^, • • • , ^^^p 1 , ordered by their energy £j(t ). 
These eigenstates evolve continuously to \ipj(t)) . The un- 
normalized adiabatic state responsible for the transport 
reads [4] 

| V'o (*) ) = cos 6 \m) A + (-l) 2 ^ sin6 \n) A 



-x£)(-l)'-|2j-2- 

3=2 



(24) 



with 



= arctan — , A = 



Note that |Vo(*o)> = \ m ) a and \^o(t m ax)) = \n) A , i.e. 
the states in which the electron is on the m th and n th 
QD, respectively. Furthermore during the entire process 
£ o{t) = holds, which ensures that no dynamic phase 
appears for the state to be transported. If 



|e (t) -£,-(*)!» 



df 



(25) 



holds and if the system at t = to is in \tpo(to)), then 
the adiabatic theorem states that the system will stay 
in \ipo(t)), provided it is a closed system. Therefore an 
electron starting in \n) A will end up in \m) A . 

To generalize this concept to the open system described 
here, we follow [25] and transform Eq. ^ with the uni- 
tary operators defined by 



uHt) \j) A 



l^i(*o)> 
\J) A 



LUL -J — 2 ' " ' ' 2 (26) 

for j < m and j > n 



to get 



7.A 



df 



dU 



Y, e J (t)\^(to))(^(t )\- l U^ — , 



~A 



N 



-Rp A + Rj2Aip A Ai 



(27) 



2 = 1 



with O 
Eq 



WOU for an operator O. If, in addition to 



( 25 1 we also assume weak coupling to the environ- 
ment [24], one can neglect the term in Eq. (27 1 as 
is shown in [35]. This is the generalization of the adia- 
batic theorem to systems described by a master equation 
of Lindblad form. Hence we have achieved the time in- 
dependence of the eigenspaces of the transformed Hamil- 
tonian H: 



dp A 



£j\ipj(to))(ipj(to)\,P 

" 2 

N 

Rp A + RY^P A ^i- 



(28) 
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For R = we get the von Neumann equation 



dp A 



E E J l^(*0)><^j(*0)] 



(29) 



and p£ (t) := (ip (to)\ P A (t) |-0p(*o)> = aHp (*o) \m) A 
is therefore a constant of motion, which ensures perfect 
transport of the electron. At t — t max we use U (t max \ to 



transform back to the Schrodinger picture. From Eq. ( 26 ) 
we find again that an electron initially in state \m) A = 



0) 



at the end of 



|-0o(*o)) will be in state \i/jo(tr, 
step two. 

For R ^ we can use Eq. (28 1 to calculate the loss 
from perfect transport 



N 



-Rp A (t) + {Mt)\A jP A (t)A J |Vo(f)> 



N 



= -R~p£S) + 



Rj 

j — l i,i'—m 



a 

TiA 



(Mt)\i) A <j'\Mt)) pfAt) (30) 



-i'j 



As we argued before, without measurements p it) is a 
constant in time in the adiabatic approximation. To first 
order in the measurements [25] we can therefore sub- 
stitute p A (t) -> p A (t ) = |^o(*o))(^o(*o)| and p A (t) 
|^o(i))(V'o(*)| ° n the right side to find 
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-R\l 
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\(M*)\- 



(31) 



where \ipo{t)) can be substituted from Eq. (24 1 to obtain 
numerical values. 

One might also ask what happens to coherences 
PfeoW := A(k\p A (t) H/taCto)) (for any fc < m or fc > n) 
during the transport. This question arises if one is in- 
terested in the transport of one part of a super position 
state, e.g. ^ (\k) A + \m) A ) -> + \n) A ). In the 

same manner as Eq. (30) and Eq. (31) we arrive at 



P&(*) 



dt 
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FIG. 3: Transfer loss during the transport of an electron along 
three (black), five (blue), seven (red), nine (green), and eleven 
(purple) QDs as a function of §, i.e. the localness of the 
measurements. The crosses show the solution from the exact 
Eq. (JsJ> for local measurments and the circles show the em- 
pirical cross-over between local and non-local measurements 
(see text). Pump and Stokes pulses are as in Eq. (371 (see 



also Fig[4][a)) with T chosen such that the additional loss due 
to non-adiabaticity is 24 x 1CP 6 for all curves, to allow rea- 
sonable comparison. This is T = 150, 196, 225, 242, 249 
for three, five, seven, nine, and eleven QDs, respectively. The 
intermediate couplings are equal one. System parameters are 
R = N, a = 0.04a. 



where again Eq. (24 1 can be substituted 



The terms in the parenthesis of Eq. (31 1 and Eq. (33 1 



become small when the measurements arc sufficiently 
non-local that they can not distinguish well between the 
QDs involved in the transport (see also Fig|3|. Since the 
loss of information increases linearly in the time t max re- 
quired for transport, it is crucial that the couplings £lj 
between the dots are as large as experimentally possible 
to give a big energy splitting which in turn allows a fast 
transport (see Eq. (25)). 

Also note that in the case of local measurements 



Eq. (132] reduces to ^iP~ A 



ence rate T2 = |r 
one obtained in Eq 



-T 2 p A where the decoher- 



2 + a — 2\f\ + a) is the same as the 
"12] ) without transport. Therefore, 
the decoherence rate of a charge qubit on a quantum dot 
rail is the same during storage as during transport by 
CTAP if subject to local measurements. 

The probability of not finding the electron in the 
desired state I 



[ after the transport, calculated from 
Tk i Eq. (31 1 is shown in Fig [3] as a function of the parameter 

2 for transport along 3, 5, 7, 9, and 11 dots. Remember 
-a(*| V'o(i)) Pki(t) (32) that 2 is the distance of the QPCs to their nearest QD 
' compared to the distance between neighboring QDs, and 

n therefore is a measure of the non-localness of the mea- 

surements performed by the QPCs. For any given a this 
parameter can be changed with d and we fix the sensitiv- 



A! 



)ity A (i\ ni \i) A of the measurements by setting a — 0.04a. 
(33) Also the rail of dots and measurement apparatuses ex- 
tends on both sides (N ^> m — n). 
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One can distinguish two regimes in Fig [3] First, for 
a > n-m (j-[g\^ f circles) the transport fidelity in- 
creases with the non-localness of the measurements. This 
is easily understood since sufficiently non-local measure- 
ments can not distinguish between the QDs involved in 
the transport. 



the Hamiltonian for N = 3 explicitly 



Second, for 



< 



(left of the circles) we find the H 



d 4 

reverse and the transport fidelity increases with the lo- 
calness of the measurements (except for transport along 
three QDs). This needs some explanation. In this regime 
the measurements are local enough to distinguish well be- 



tween the states \m) A and 



(during CTAP most of 



the population is found on these two states) and hence 
one can not expect a further decrease of the fidelity with 
the localness of the measurements. To account for the 
apparent increase of the fidelity with localness we need 
another argument. For very local measurements, sites n 
and m almost exclusively contribute to the decoherence 
whereas if the measurements are slightly non-local, sites 
in the near neighborhood of n and m also contribute to 
decoherence, therefore increasing the over all decoher- 
ence. 

The small decrease of the fidelity for extremely local 
measurements might be due to the small populations on 
QDs other than m and n. 

It is also evident from Fig |3j that if local dephasing is 
the main source of decoherence, it is best to transport 
long distances in one step as is seen e.g. for the transfer 
1 1) A — > 1 11)^. By doing so the chance of not finding the 
electron on the desired state after the transport is about 
11 * 1CT 3 (see left end of Fig [3]). This number increases 
to about 5* (4* 10 -3 ) when we first transfer to |3) A , then 
to |5) A , \1) A , |9) A , and finally to |11) A . On the contrary, 
if dephasing tends to be more global, a better transfer 
rate is achieved by breaking a long distance into several 
smaller ones. 



B. Coupling to TLSs 



We now apply CTAP to QDs coupled to TLSs as shown 
in Fig [2] Hence the driving field Ha from Eq. ( 23 ) has 
to be added to H int from Eq. ( |13| |, responsible for the 
QD-TLS-coupling. We first neglect the coupling to the 
measurement apparatus and point out it's inclusion at 
the end of the subsection. With this all couplings behave 
local and we can restrict ourself to study only the QDs of 
the rail which are involved in CTAP, i.e. N = n — m + 1. 

To understand how CTAP works for this system, we 
note that the Hamiltonian H — Hi nt + Ha can be writ- 
ten in a block diagonal form, with 2 N blocks, each block 
having dimension N [26] . Each block gives the Hamil- 
tonian of the electron on the rail of QDs conditioned on 
the state of the TLSs. As an example we show part of 
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where the first block is for all TLSs in the |0) state, the 
second for the first and second TLSs in |0) and the third 
TLS in |1), and so on. 

We can now study each block individually because 
there is no coupling between them. If we want to trans- 
port the electron from the first QD to the third QD, re- 
gardless in which state or superposition the TLSs are, we 
have to make sure that CTAP works in each block. That 
is the state \1) A |0) |0) |0) should adiabatically evolve to 
|3) A |0) |0) |0) as well as \1) A |0) |0) |1) should evolve to 
|3) A |0) |0) |1) and so on. If this is achieved, then an 
electron initially found in 11)^ will evolve in the pure 
state 1 3)^. During the transition the reduced state 
p A = r FrTLSs[p\ is not only in a superposition of 11)^, 
1 2) j_and 13)^, but in a real mixture of these states (see 
Fig[5|a)). This is because depending on the state of the 
TLSs, the electron might start moving earlier or later, 
which generally results in cntaglcmcnt between the rail of 
QDs and the TLSs. However, once CTAP is completed, 
the electron will be found in the desired state 13)^. 

Hence, for CTAP to be applicable despite coupling to 
TLSs, we only have to make sure that it works for a 
Hamiltonian of the form (where we return to an arbitrary 
long rail of dots) 



H = 



( ±Xm ftm 
ftm ftm+1 
ftjn+1 ±Xm+ 2 





V 



±x, 



\ 




ftn-1 
±Xn / 



(34) 



instead of the simpler form Eq. ( 23 1 . CTAP should work 
for all possible permutations of + and — , since each per- 
mutation represents one block of the Hamiltonian corre- 
sponding to a certain combination of environmental TLS 
states. The crucial requirement for this is again Eq. (25 1, 



but this time it can not be fulfilled by just performing 
the transport sufficiently slowly. One also has to make 
sure, that no level crossings occur, i.e. Eo(t) — Ej(t) ^ 
during step two. As before, £o(^) is the energy of the 
state \tj}o(t)} which we want to evolve adiabatically from 
IV'o(to)) = \m) A to \i/j (t max )) = \n) A . Note that un- 
like in the previous subsection, £o(t) ^ and depends on 



the state of the TLSs. The situation becomes even more 
complicated when anti-crossings are considered which ap- 
pear quite naturally (see Fig [3]) . At an anti crossing two 
energy levels get so close that adiabaticity can not be 
achieved with any reasonable transfer times. 

Furthermore it might happen that the adiabatic state 
\ipo(t)) does not connect \m) A to \n) A , but instead to an- 
other state \ipo(t m ax)} = \j) a w hh j ^ n. An example of 
this is shown in Fig [4] (blue curves) where the population 
returns to the original QD. 

To calculate the energies is generally only possible nu- 
merically and an exact treatment when level crossings or 
anti-crossings appear is highly non-trivial. Some work 
on this issue is done in [T3], but the possibility of anti- 
crossings is not considered there (see also [H]). However, 
some qualitative statements can be said. If the energy 
level e (t) of the adiabatic state is at all times during 
step two in the center of all energy levels involved in the 
transport, i.e. 



£ (t)>E j (t) for j 



m — n 



e (t) <£,•(«) for j = l,...,— — , (35) 



then we can be sure that no level crossings occur. In this 
case the process during step two is qualitatively the same 
as without coupling to the TLSs. Since condition (35 1 is 
valid with \j = (previous subsection), it can always be 



achieved with sufficiently high driving fields fij, as then 
the comparatively small couplings \j do not influence the 
qualitative eigenvalue structure. 

Lets apply this condition to the block of the Hamilto- 
nian in which we have minus for Xm and pluses for all 
other coupling constants. Clearly we need some strong 
couplings flj(t a ) to ensure that £0(^0) — ~Xm is n °t 
the lowest energy, but the centrefold one (see Fig Qc) 
where this crossing occurs at i ~ 5 before the trans- 
port starts at to w 20). For n — m = 2 we find 
^n-i(*o) > (Xi - Xa){Xi ~ X2)- For transport along 
many QDs (n — m 3> 2) we get the approximate in- 
equality f2j(fo) ~ cmax(xt)(fi- m) for some constant c, 
which states that for given maximal Qj (due to the ex- 
perimental setup) and given \j 1 there exists a maximum 
number n — m for which we can use CTAP. For transport 
along more QDs one has to break up the transport into 
several smaller sections and apply CTAP successively to 
each section. The shorter the rail gets (compared to the 
maximal length) , the further we can stay away from any 
level crossing and the faster we can transport the electron 
without violating the adiabaticity condition (251. 



As example we show in Fig [4] the cigcncncrgies and 
populations as a function of time for a system with five 
QDs and constant Xi — X f° r different values of \. We 
consider the block of the Hamiltonian acting in the sub- 
space where the first TLS is in the state |0) and the other 
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FIG. 4: (a): Coupling between QDs (dashed: fii, dotted: ^2,3, 
solid: ^4), (b)-(e): energy levels of Hamiltonian Eq. (361 and 



X = 0, 0.15, 0.3, 0.45 (see text), and (f)-(j): populations 
of |1) 4 - |5} A obtained by the numerical integration of the 
Schrodinger equation. The colors black, red, green and blue 
correspond to x = 0, 0.15, 0.3, 0.45, respectively. What looks 
like level crossings in (c) and (d) are actually anti-crossings. 
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four TLSs are in the state II), which reads 



Ha 



( -x \ 

ftp + x n 2 o o » 

o o 2 +x n 3 o 

o o q 3 + x n s 

V o o o n s + X J 



(36) 



Of course, if we assume the TLSs to be in a statisti- 
cal mixture one would have to calculate the populations 
for all blocks of the Hamiltonian and average over them. 
Since all qualitative features can be seen with this par- 
ticular block we restrict the discussion to this one for 
pedagogic reasons. We use 2 — ^3 = 1 and Gaussian 
pump and Stokes pulses 



(t - 3T/4) 2 



ftp = exp ^ — 
fls — exp 



(T/4)2 
(t - T/4) 2 



{T/Af 



(37) 



with the parameter T = 150. Because of the Gaussian 
pulse shapes the distinction of steps one, two and three is 
not sharp, any more. Looking at Fig[4](a) we can approx- 
imately say that step one is for —50 < t < 20, step two is 
for 20 < t < 130 and step three is for 130 < t < 200. T is 
chosen long enough to achieve very good transfer fidelity 
(0.9996) for x = 0, i.e. without coupling to TLSs. For 
X = 0.15 there appears a level anti-crossing at which the 
energy gap (0.0002) is so small, that in can be treated as 
a crossing. However, this happens during step one before 
the actual population transfer starts and does therefore 
not disturb adiabaticity much. The smaller population 
transfer (0.975) is because during step two the energy 
level spacing is decreased compared to the previous case. 
Therefore, the transfer rate can be increased arbitrarily 
close to one by choosing larger T [27). Enlarging % = 0.3 
moves the anti-crossing to step two and adiabaticity can 
not be achieved any more. The result is wild oscilla- 
tions of all populations and a low transfer rate. Finally 
at x = 0-45 the anti-crossing disappears and adiabatic- 
ity is restored which is evident in the lack of oscillations 
in the populations. But because the eigenenergy -0.15 
at t = — 00 connects to the same energy at t = 00, the 
population returns to the original state \m) A . 

This result should be quite surprising. Up to relatively 
large couplings to TLSs, — < 0.2 (for m—n = 4) adi- 
abatic transfer can be achieved and the transfer fidelity 
very quickly approaches one with sufficiently large trans- 
fer times. This is contrary to the Markovian dephasing 
studied in the previous subsection, where the transfer 
loss increases with transfer time once the transfer time is 
long enough to ensure adiabaticity. 

One disturbing effect introduced by the coupling to 
TLSs is that £q ^ and hence we get a dynamical phase 



If : 



dteo(t) 



(38) 
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FIG. 5: The purity of the reduced density matrix p A (t) dur- 
ing the transfer 1) 4 — > |3) A . 

(a) : Transport of a position eigenstate. \ip(t = 0)) = |1) A and 
Xi = 0.02 (solid), 0.05 (dashed), 0.15 (dotted). Because of the 
dependance of \ipo(t)) on the state of the TLSs, the electron 
gets entangled with the TLSs during transport. However, this 
entanglement disappears with the completion of CTAP and 
the electron on the quantum dot rail finishes the transport in 
the desired pure state. 

(b) : Transport of half of a superposition state. \tf)(t = 0)) = 
i (\0} A + \1) A ) and x* = 0.005 (solid), 0.02 (dashed), 
0.1 (dotted). The transported half of the superposition state 
picks up a dynamical phase which depends on the state of the 
TLSs. Hence the electron gets entangled with the TLSs and 
finishes CTAP in a statistical mixture. 



and worse even, this phase depends highly on the state 
of the TLSs. If the electron starts in a superposition 
^(\k) A + \m) A ) it will not end up in ^ (\k) A + \n) A ) 
but in a real mixture of \k) A and \n) A . That is, dur- 
ing the transport, the electron gets entangeled with the 
TLSs. This should not be surprising since we saw a simi- 
lar behavior already in subsection |II B| without transport 
(see Eq. (21 1). However, CTAP prevents the lost infor- 



mation from returning back to the electron periodically. 
This can be seen in Fig [5jb) for transport along three 
QDs, Xt = X, and fl P , fl s as in Eq. ^ with T = 150. 
The initial states Pj(t ) = 1 1 of the TLSs are taken to 
be a complete mixture. As a measure of purity we use 



Tr 



On the other hand, as described earlier this 



section, if the electron starts in 11)^ it will get transfered 
to the pure state (3)^ despite being in a real mixture 
during the transfer as is shown in Fig[5ja). 

Now we briefly discuss the inclusion of measurements 
once again. In the previous subsection the adiabatic the- 
orem was assumed to hold and therefore we first have to 
make sure that it holds also with coupling to TLSs. This 
is best done by solving the Schrodiner equation without 
measurements as was done to get Fig [4] The solution 
shows whether the transformation is adiabatic (almost 
no fast oscillations) or not (much oscillations). If not, 
the driving fields f2 3 as well as the transfer time can be 
increased, until the solution shows adiabatic behavior. 
Then equations 



-(33) derived in the previous subsec- 
tion can be applied here as well, if one uses the appro- 
priate eigenstate |V'o(^)) which depends on the state of 
the TLSs. However, since \ipo(t)) can only be calculated 
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numerically and the effects of dephasing do not depend 
much on the exact form of this state, one might prefer 



to use the unperturbed eigenstate from Eq. (24 1 as in 



subsection |III A| It is therefore justified to examine deco- 
herence effects arising due to measurements and due to 
coupling to TLSs separately. 



IV. DISCUSSION AND CONCLUSION 

While storing information as a charge qubit on a quan- 
tum dot rail, some of it will leak because of decoherence 
due to measurements. Information loss from decoherence 
arising from coupling to TLSs, on the contrary, returns 
periodically. 

On the other hand, while transporting information 
along a rail of QDs with CTAP, information lost to TLSs 
will not return to the electron on the QDs. However, only 
information of the phase will be lost since in the adia- 
batic limit one can achieve perfect population transport 
despite coupling to TLSs, provided it is not too strong 
compared to the driving field flj and the QDs involved 
are few enough. This is different for decoherence through 
measurements when there is always a chance of loosing 
the electron during the transport, which can be mini- 
mized by doing the transport fast. 

Summarizing we can say, if we can reduce Markovian 
and non-Markovian dephasing sufficiently to be able to 
store a charge qubit on a rail of QDs, then we can also 
transport it using CTAP. 

It is often claimed that CTAP is relatively robust to 
decoherence because the QDs between the initial and the 



final QDs are barely populated. This is in contrast to our 
results because of the following reasoning. To achieve a 
small population on the intermediate QDs the coupling 
between these has to be large compared to the Stokes 
and pump pulses, as seen in Eq. p4| . On the other hand 
we find that all couplings should be as large as possi- 
ble (for both, Markovian and non-Markovian decoher- 
ence) to achieve big energy spacings which in turn allow 
fast transport. Therefore it is reasonable to choose large 
Stokes and pump pulses and allow some population on 
the intermediate QDs. The above statement might stem 
from analyses of STIRAP where the intermediate atomic 
states are exposed to decay, but this is not the case for 
CTAP. 

Although the non-Markovian noise studied here is 
fairly simple and might not be very realistic in an ex- 
periment, the striking differences we find compared to 
Markovian dephasing indicates that it might be worth 
studying more specific examples of this class of decoher- 
ence. 

Another approach to quantum computing is the use of 
the spin degree of freedom of an electron as a qubit. In 
this case the differences between the influence of Marko- 
vian and non-Markovian dephasing on transport with 
CTAP might become even more apparent. Because the 
coupling of the spin degree to the environment is ex- 
pected to be much weaker than the coupling of the 
charge, the dynamic phase acquired during the transport 
due to coupling to TLS's should be fairly independent of 
the spin and therefore an undetectable global phase even 
for superposition of spin states. 
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